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In this note we mainly explore the possibility of measuring the action of the intrinsic 
gravitomagnetic field of the rotating Earth on the orbital motion of the Moon with the 
Lunar Laser Ranging (LLR) technique. Expected improvements in it should push the 
precision in measuring the Earth-Moon range to the mm level; the present-day Root- 
Mean-Square (RMS) accuracy in reconstructing the radial component of the lunar orbit 
is about 2 cm; its harmonic terms can be determined at the mm level. The current 
uncertainty in measuring the lunar precession rates is about 10~^ milliarcseconds per 
year. The Lense-Thirring secular, i.e. averaged over one orbital period, precessions of 
the node and the perigee of the Moon induced by the Earth's spin angular momentum 
amount to 10~^ milliarcseconds per year yielding transverse and normal shifts of 10~^ — 
10~^ cm yr~^. In the radial direction there is only a short-period, i.e. non-averaged over 
one orbital revolution, oscillation with an amplitude of 10~^ m. Major limitations come 
also from some systematic errors induced by orbital perturbations of classical origin like, 
e.g., the secular precessions induced by the Sun and the oblateness of the Moon whose 
mismodelled parts are several times larger than the Lense-Thirring signal. The present 
analysis holds also for the Lue-Starkman perigee precession due to the multidimensional 
braneworld model by Dvali, Gabadadze and Porrati (DGP); indeed, it amounts to about 
5 X 10""' milliarcseconds per year. 
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1. Introduction 

In the framework of the linearized weak-fiel d an d slow-motion approximation of 
general relativity, the gravitomagnetic effectJiEl ^j-g in duc ed by the off-diagonal 
components goi, i = 1, 2, 3 of the space-time metric tensoJ^^ which are proportional 
to the components of the matter current density of the source ii = pvi . 

There are essentially two types of mass currents in gravity^. The first type is 
induced by the rotation of the matter source around its center of mass and generates 
the intrinsic gravitomagnetic field which is closely related to the proper angular mo- 
mentum S (i.e. spin) of the rotating body. The other type is due to the translational 
motion of the source and is responsible for the extrinsic gravitomagnetic field. 

A debate has recently arisen concerning the possibility of measuring some ex- 
trinsic gravitomagnetic orbital effects affecting the motion of the Earth-Moon sys- 
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tern in the Sun's field with the Lunar Laser Ranging (LLR) 

techniqujsmsmoin]. 

Another test of extrinsic gravitomagnetism concerning the deflection of electromag- 
netic waves by Jupiter in its o rbit al motion has been performed in a dedicated 
radio-interferometric experiment^]. 

In this brief note we wish to consider in some details the possibility of measuring 
with LLR an effect induced by the intrinsic gravitomagnetic field of the spinning 
Earth 

G\3r{r ■ S) -r^S] 

through the non-central, Lorcntz-like acceleration 

a = -2(^)xB, (2) 

on the orbital motion of the Moon around the Earth. In eq. ([T|) and eq. G is 
the Newtonian gravitational constant, c is the speed of light in vacuum, S is the 
Earth's spin angular momentum and v is the velocity of the Moon with respect to 
the Earth. The orbital feature we are interested in consists of the Lense-Thirring 
precessions of the longitude of the ascending node and the argument of pericentre 
w of the orbit of a test particljl^l Twenty years ago, Bertotti in Ref. [HI wrote 
that the Lense-Thirring effect for the Moon was, at that time, too small to be 
detected; according to CiufolinP^, intrinsic gravitomagnetism is still unmeasurable 
with the lunar orbit. Instead, the possibility of measuring it in view of the expected 
forthcoming improvements in LLR has recently been envisaged by Miiller et al. 
in Ref. [15] and in Ref. [161 Kopeikin et al. in Ref. [TTj retain highly plausible its 
measurement. An overview of other attempts to measure the Lense-Thirring effect 
in various Solar System scenarios with natural and artificial test particles can be 
found in Ref.fTSl In particular, for the LAGEOS-LAGEOS II test in the gravitational 
field of the Earth see Ref. \19\ and Ref. 20l for the perspectives on measuring the 
solar intrinsic gravitomagentic field with the inner planets of the Solar System see 
Ref. [2TJ for the test performed with the Mars Global Surveyor in the field of Mars 
see Ref. [22] and Ref. [23l Another effect induced by the in trinsic gravitomagentic 
field of the Earth is the precession of orbiting gyroscopeJ^lES] currently under 
measurement by the GP-B mission; see on the WEB http:/ /einstein. Stanford. edil/l 
Ref. [26] and Ref. [27] 

It is interesting to note that our analysis is equally valid also for the anoma- 

952 

lous Lue-Starkmart^ perigee precession predicted in the framework of the multi- 
dimensional braneworld model of modified gravity put forth by Dvali, Gabadadze 
and Porrat^22] (DGP) to explain the observed acceleration of the Universe without 
resorting to dark energy; indeed, as we will see, the magnitude of such an effect is 

oo on oi 

the same as the Lense-Thirring one for the Moon. Several researchers ° •J" •J^ ar- 
gued that it might be possible to measure the Lue-Starkman precession with LLR 
in view of the expected improvements in such a technique. 
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The physical and geocentric orbital parameters of the Moon are listed in In 
Table m 



Table 1. Physical and geocentric orbital parameters of the MoorJSSESl. 
The gravity field adopted is the LP150Q solution (See on the WEB 
|http://pd s-geosciences.wustl.edu/geo/lp-l-rss-5-gravity-vl/lp_1001/sha/jgll50ql.lbl^. 



Parameter 


Value 


Units 


m mass 


7.349 X 10^2 


kg 


S proper angular momentum 


2.32 X 10^9 


kg m^s~^ 


Gm 


4.902801076x1012 


]jj3 g-2 


R radius 


1.738 X 10'5 


m 


a proper angular velocity 


2.66x10-6 


rad s— 


normalized moment of inertia 
J2 mass quadrupole moment 


0.3932 




2.0326 X 10"* 




SJ2 uncertainty in the mass quadrupole moment 


1 X 10"** 




a semi-major axis 


3.84400 X 10** 


m 


I mean inclination to the Earth's equator 


23.5 


deg 


e eccentricity 


0.0549 





2. The Lense-Thirring effect on the lunar orbit 

By assuming a suitably constructed geocentric equatorial frame, it turns out that the 
node and the perigee of the Moon undergo the Lense-Thirring secular precessions 

' ^^LT ^ c2a3^fff2)3/2 = 0.001 mas yr-\ 

(3) 

QGSe^cosI nno —1 

^ WLT = - e2a3(i_e'^)3/2 = ^0.003 mas yr \ 

where a, e, / are the semi-major axis, the eccentricity and the inclination to the 
Earth's equator of the Moon's orbit: ruas yr"^ stands for milliarcseconds per year. 
We used — 5.85 x 10'^^ kg m^ s~^'311_ xhe Lue-StarkmaiJ^ pericentre precession 
is just aboulE =f0.005 mas yr"^. 

Since the ratio of the mass of the Moon to that of the Earth 
0.0123000383, one may argue that eq. which has been derived for a test-particle 
like, e.g., an artificial satellite, does not apply to the Earth- Moon system. The in- 
trinsic gravitomagnetic spin-orbit effects in the case of a two-body system with 
arbitrary masses ttia and tob and spins and have been derived by Barker 
and O'Connell in Ref. 36, Damour in Ref. [37l Wex in Ref. [38l for the sake of sim- 
plicity, we will reason in terms of the node. In this case, the total node precession 



^The minus sign is related to the standard Friedmann-Lemaitre-Robertson- Walker (FLRW) 
branch, while the plus sign is related to the self-accelerated branch which s hou ld be able to explain 
the observed acceleration of the Universe without resorting to dark energjl2Q| 
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f^tot accounts for the spin-orbit contributions of both bodies and also for a spin-spin 
term. The expression of the node precession of a body A is^*^ "^"^ 

so that 



2c2 J a'^(l — e^)'^/^ mA ' mA + TOb ' 

Jitot = i^A + f^B- (5) 



Let us pose 

mA = TOg, = Af, niB = mMoon = m; (6) 
thus, it is possible to obtain 

l^l + i^j ,2^3(1 2)3/2' (7) 

"Moon = (l + ^jj^^^; (8) 

recall that fi is the Moon/Earth mass ratio. It results that the precession of eq. 
([7]) is larger than the Lense-Thirring one of eq. ^ by the multiplicative factor 
(l + ln) = 1.0092 yielding an error of 10 mas yr , which is completely negligible 
(see Section [3]). Concerning the precession due to the lunar spin, we have 

^^Moon / 3 -|- 4^ \ 1 S'm oon 1 r» — 3 /r\\ 

— ' ' = 2 X 10 , (9) 



i.e. it is of the order of 2 x 10^® mas yr ~\ which is negligible as well. The amplitude 
of the spin-spin term is proportional tJSSETEHI 



• 3 GM{1 + ^) 1 5© S'Moon „ -,„_9 _i /,„x 

Thus, we can conclude that the Lense-Thirring approximation is fully adequate for 
the Earth-Moon system. 



3. Some sources of error 

Let us now examine some sources os systematic errors. In regard to the potentially 
corrupting action of the mismodelling in the even = 2,4,6,...) zonal (m = 0) 
harmonic coefficients Jt of the multipolar expansion of the Newtonian part of the 
Earth's gravitational potential, wh ich is not the most important source of aliasing 
precessions in the case of the MoorP2l, only (5 J2 would be of some concern. Indeed, 
the mismodelled secular precessions induced by it on the lunar node and perigee 
amount tcQ —2.67 x 10"'* mas yr~^ and 5.3 x 10"'' mas yr^^, respectively; the 
impact of the other higher degree even zonals is negligible being < 10~* mas yr~'. 

''The calibrated errors 5Ji of the EIGEN-CGOIC Earth gravity field solutioiiSSl were used. 
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As in the case of the spins, also the asphericity of the Moon has to be taken into 
account) according to 



n 



Moon 



3 TT-Moon 

2 



cosFJ2^°°" /i?, 



(1 



Moon 



(11) 



where riMoon = \/ GM(1 + ij)/a^ is the lunar mean motion and F is the angle 
between the orbital angular momentum and the Moon's spin angular momentum 
S^Moon; it is about 3.61 deg since the spin axis of the Moon is tilted by 1.54 de g to 
the ecliptic and the orbital plane has an inclination of 5.15 deg to the ecliptic^. 
Table [Hand eq. (fTTj) yield a mismodelled node precession due to i5J|*°™ of about 
0.006 mas yr~^, which is 6 times larger than the Lense-Thirring rate. For other 
sources of systematic errors induced by gravitational and even non-gravitational^ 
perturbations see Ref. ^2^ and references therein, especially Ref. |41] Among the N- 
body gravitational perturbations, the largest ones are due to the Sun's attraction. 
In order to get an order-of-magnitude evaluation of their mismodelling, let us note 
that some of such effects are proportional to n^/riMoon; e.g. the node rate, referred 
to the equator, iJ 



SGOJIqcos/ /3 

4amnMoon 



-5 X 10^ mas yr ^, 



(12) 



where e = 23.439 deg is the obliquity of the ecliptic. SincJ^ SGdyt(. 



.cermng the 
j44l45ll5ll6l 



8 X 10^ m'^ s~2, we can assume a bias of « 0.07 mas yr~^ which 
is 70 times larger than the Lense-Thirring precession. 

Let us, now, consider the precision of LLR in reconstructing the lunar orbit 
with respect to the Lense-Thirring effect. Concerning__the^ precision in measuring 

0.1 mas yr~^, i.e. it is 



the lunar precession rates, it amounts to abou 
two orders of magnitude larger than the Lense-Thirring precessions of eq. The 
orbital perturbations experienced by a test particle are usually decomposed along 
three orthogonal directions of a frame co-moving with it; they are named radial 
R (along the radius vector), transverse T (orthogonal to the radius vector, in the 
osculating orbital plane) and normal N (along the orbital angular momentum, out of 
the osculating orbital plane). According to, e.g., Ref. HBl the R—T—N perturbations 
can be expressed in terms of the shifts in the Keplerian orbital elements as 



Ai? 



AT = 



AiV 



(Aa) 



2 ■ 



AX + Acj + cos/Ar2 + y^Kef+JeAMf 



(13) 



2 



(sin mny 
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where M is the mean anomaly. The lunar Lense-Thirring shifts after one year are, 
thus 

' Ai?LT = 0, 

< ATlt = a\/l + T (^^LT + cos/A^2lt) = -0.38 cm, (14) 



^ AiVLT = ay 1 + T sin/Ar^LT = 0.07 cm. 

It is important to note that there is no Lense-Thirring secular signature in the 
Earth-Moon radial motion on which all of the efforts of LLR community have been 
concentrated so far. It can be shown that a short-period, i.e. not averaged over one 
orbital revolution, radial signal exists; it is proportional to 

Ar oc = 2 X 10-5 m, (15) 

c na 

which is too small to be detected since the present-day accuracy in estimating the 
amplitudes of radial harmonic signals is of the order of mrrP. Major limitations 
come from the post-fit Root-Mean-Square (RMS) accuracy with which the lunar 
orbit can be reconstructed; the present-day accuracy is about 2 cm in the radial 
direction R along the centers-of-mass of the Earth and the Moorff^. Improvements 
in the precision of the Earth-Moon ran ging of the order of 1 mm are expected in the 
near future with the APOLLO progranP^EH Recently, sub-centimeter precision in 
determining range distances between a lase r on the Earth's surface and a retro- 
reflector on the Moon has been achievecf^. However, it must be considered that 
the RMS accuracy in the T and N directions is likely worse than in R. 



4. Conclusions 

In this note we have examined the possibility of measuring the action of the intrinsic 
gravitomagnetic field of the spinning Earth on the lunar orbital motion with the 
LLR technique. After showing that the Lense-Thirring approximation is adequate 
for the Earth-Moon system, we found that the Lense-Thirring secular precessions of 
the Moon's node and the perigee induced by the Earth's spin angular momentum 
are of the order of 10"'^ mas yr~^ corresponding to transverse and normal secular 
shifts of 10^^ — IQ-^ cm yr"^. The intrinsic gravitomagnetic field of the Earth 
does not secularly affect the radial component of the Moon's orbit; a short-period, 
i.e. not averaged over one orbital revolution, radial oscillation is present, but its 
amplitude is of the order of 10"^ m. The current RMS accuracy in reconstructing 
the lunar orbit is of the order of cm in the radial direction; the harmonic components 
can be determined at the mm level. Forthcoming expected improvements in LLR 
should allow to reach the mm precision in the Earth- Moon ranging. The present-day 
accuracy in measuring the lunar precessional rate is of the order of lO^^ mas yr^^. 
Major limitations come also from some orbital perturbations of classical origin like. 
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e.g., the secular node precessions induced by the Sun and the oblateness of the 
Moon which act as systematic errors and whose mismodelled parts are up to 70 
times larger than the Lense-Thirring effects. As a consequence of our analysis, we 
are skeptical concerning the possibility of measuring intrinsic gravitomagnetism with 
LLR in a foreseeable future. The same conclusion holds also for the Lue-Starkman 
perigee precession predicted in the framework of the multidimensional braneworld 
DGP model of modified gravity; indeed, it is as large as the Lense-Thirring one for 
the Moon. 
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